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We prove that the threshold detection efficiency for a loophole-free Bell experiment using an 
n-qubit Greenberger-Horne-Zeilinger state and the correlations appearing in the n-partite Mermin 
inequality is n/(2n — 2). If the detection efficiency is equal to or lower than this value, there are local 
hidden variable models that can simulate all the quantum predictions. If the detection efficiency 
is above this value, there is no local hidden variable model that can simulate all the quantum 
predictions. 
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Quantum nonlocality is the impossibility of reproduc- 
ing the quantum correlations between the results of dis- 
tant measurements using local hidden variable (LHV) 
theories. This impossibility is shown either by the vi- 
olation of a Bell inequality [T] or the impossibility of as- 
cribing predefined results which simultaneously satisfy 
several predictions of quantum mechanics [^. Although 
quantum nonlocality is intimately related to entangle- 
ment [3] , security of quantum cryptography J? , and com- 
munication complexity ^ , there is as of yet no loophole- 
free quantum nonlocality experiment. A particulary im- 
portant problem is the detection loophole. It occurs when 
the imperfect efficiency of the detectors leaves room for 
LHV theories in which undetected events can occur due 
to local hidden instructions rather than to imperfections 
[B]. An appropriate measure of the quantum nonlocality 
of a given quantum state and Bell inequality is, therefore, 
the minimum detection efficiency required for a loophole- 
free Bell experiment, ?7crit- It is defined as the value of the 
ratio between detected and emitted particles such that, if 
V ^ ^ciiti there is a LHV theory reproducing the predic- 
tions of quantum mechanics, but no such LHV theories 
exist if ?7 > r^crit • The value of r]„it is known for some sce- 
narios [3 m El Uni HH [H] , and some general bounds have 
been obtained [13]. Curiously, r^crit was still unknown for 
a very important scenario. 

Eighteen years ago, Mermin and others discovered 
the first example of a Bell inequality with a violation 
that grows exponentially with the number n of parti- 
cles [21 [TS]. Specifically, they show that the n-qubit 
Greenberger-Horne-Zeilinger (GHZ) state |GHZ„) ^ vi- 
olates a n-partite Bell inequality by an amount that 
grows as 2("-i)/2 [14l[15]. If instead of a pure |GHZ„) we 
have a noisy one, 1^|GHZ„)(GHZ„| + {1 - V)t/T\ then 
the minimum value of V required to observe a violation 
is Vcrit — 2^^^")/^. Later, Werner and Wolf proved that 
the Mermin inequality is the two-setting correlation Bell 
inequality "which can be violated by the widest margin 
in quantum theory . . . [and] is the only one for which the 
maximal violation 2*^""^)/^ is attained" [TB] . 



There have been several attempts to obtain rycrit for 
the Mermin inequality: Braunstein and Mann showed 
that 7?crit > 2(i-")/2" for n odd, and r/crit > 2(2-")/2" 
for n even ^7|; Brassard, Broadbent, and Tapp showed 
that rycrit < 2*^^~"^/" [18j; and Larsson proved that 
'7crit = 3/4 for n = 3 0. However, no formula was 
known for arbitrary n |19j . In this Letter we prove that 
'7crit = n/{2n — 2). In addition, we obtain numerically 
the relation between ?7crit and Vcxit for several values of 
n. 

The Mermin inequality is based on the GHZ proof of 
Bell's theorem [2" . It shows the impossibility of assigning 
predefined values —1 or 1 to local observables, simulta- 
neously satisfying several perfect correlations predicted 
by quantum mechanics. The scenario for the GHZ proof 
is the following. A system composed of n > 3 particles 
is initially prepared in the state |GHZ„). Each particle 
moves away to a distant space-time region where an ob- 
server measures randomly either Xi or Zi, where X and 
Z denote the Pauli matrices <Tx and ct^ , and i denotes par- 
ticle i. Local measurements on particle i are assumed to 
be spacelike separated from the choices of measurements 
made on all other particles. 

The n-qubit GHZ state is the unique simultaneous 
eigenstate that satisfies 

5,|GHZ„) = |GHZ„), fori-l,...,n, (1) 

where 

n 

g, = X,I^Z, (2) 

are the generators of the stabilizer group of the GHZ 
state, defined as the set {sj}^^^ of all products of the 
generators. The perfect correlations of the GHZ state 
are 

(GHZ„|s,|GHZ„) = 1, for i = 1,...,2". (3) 

For n odd, the GHZ proof is as follows. Any LHV 
theory assigning predefined values —1 or 1 to and Zi 
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in agreement with the quantum predictions given by ([s]) 
must satisfy 



1, for j — 1, . 



(4) 



However, if we take into consideration only the 2"^^ pre- 
dictions Q involving stabilizing operators sj that are 
products of an odd number of generators, and assume 
predefined values either —1 or 1, then it so happens that, 
at most, only 2""^ + 2("-3)/2 out of these 2""! pre- 
dictions are satisfied. For the remaining predictions of 
quantum mechanics, the corresponding prediction of the 
LHV model is the opposite (i.e., Sj = — 1); the reason for 
this behavior will be explained below [see (ii)]. There- 
fore, this discrepancy between quantum mechanics and 
LHV theories can be reformulated as a violation of a Bell 
inequality. Any LHV theory must satisfy the following 
inequality: 



|/3„| < 2("-i)/2^ 



where the Bell operator 



/3« = 



(5) 



(6) 



is the sum of all stabilizing operators which are products 
of an odd number of generators. Inequality (|5| is the 
Mermin inequality for n odd |14j . On the other hand, 
the n-qubit GHZ state satisfies 



(GHZ„|/3„|GHZ„) =2"- 



(7) 



and therefore violates the Mermin inequality ([5| by an 
amount that grows as 2'"~^)/^ [M]. For n even, the Mer- 
min inequality is not violated by 2^""^)/^. The equivalent 
inequality was found by Ardehali |15j . For an explana- 
tion of the Ardehali inequality in terms of stabilizers of 
the GHZ state, see [20]. For simplicity's sake, we will 
focus on the Mermin inequality for n odd. Our proof 
works similarly for n even when we consider the Arde- 
hali inequality. 

Following |10j . we include the detector inefliciency in 
the LHV model, so that the model consists of a set of 
instructions telling the n particles what to do if X or 
Z are measured. For a given particle, the only possible 
instructions are "give a detection (—1 or 1)" or "do not 
give a detection." 

P{Xi) is the probability that particle i is detected (giv- 
ing either —1 or 1) when is measured. P{Xi\Xj) is 
the probability that particle i is detected when Xi is mea- 
sured if particle j ^ i is detected when Xj is measured. 
P{Xi\XjZk) is the probability that particle i is detected 
when Xi is measured if particle j (j ^ i) is detected 
when Xj is measured and particle k {i ^ k ^ j ^ i) \& 
detected when Zk is measured. Analogously, P{XiZj\Zk) 
is the probability that particle i is detected when Xi is 



measured and particle j {j ^ i) is detected when Zj is 
measured if particle k (i ^ k ^ j ^ i) is detected when 
Zk is measured. 

In our LHV models, measurement results are prede- 
fined and are independent of the measurements on other 
particles. In addition, they must satisfy some restrictions 
dictated by the expected (and testable) behavior of the 
detectors and the properties of the n-qubit GHZ state 
for the measurements involved in a test of the Mermin 
inequality. Specifically, the following assumptions lead to 
the following restrictions. 

(i) All detectors have equal, constant detection ef- 
ficiency. The efficiency is the same when X or Z 
are measured. The detection errors are independent. 
The detectors have no dark counts. — From these as- 
sumptions it follows that P{Ai) = p, VA e {X,Z} 
and Vz e {l,2,...,n}, and P{A,, . . . , Bj\Ck, . . . , Di) = 
p'', yA,...,B,C,...,D e {X,Z} and V (different) 
i, . . . , j,k, . . . ,1 £ {1,2,... n}; r is the number of ele- 



ments in Ai 



,Bj. That is, the different probabilities 



must be symmetric under particle permutation and un- 
der the permutation of Xi and Zi. If p is the minimum 
over all possible LHV models, then rycrit — P- 

(ii) Compatibility with the statistical predictions of 
quantum mechanics for the Mermin inequality using the 
n-qubit GHZ state. — Each of the terms obtained by ex- 
panding ([6]), e.g., XiZ2Z^ . . . Zn, represents an experi- 
mental configuration required to test inequality (|5]). We 
will also consider configurations obtained for the previous 
ones by selecting measurements on subsets containing an 
odd number of particles, e.g., XiZ2Z^. According to the 
predictions of quantum mechanics for the n-qubit GHZ 
state, in each of these experimental configurations, when 
an odd number 3 < g < n of particles are detected, the 
corresponding results must satisfy 

XiZ^Zk ■ ■ ■ Z„ = ZiXjZk ■ ■ ■ Z„ = ZiZ^Xk ■ ■ ■ z„ = ■ ■ ■ 

(8) 



— ZiZiZi. ■ ■ ■ X„ — —XiXiXh ■ ■ ■ x„. 



In addition, if g ^ n, then (|8| must equal Z^j^-y . . . Z„. 
Therefore, depending on the result, —1 or 1, of the prod- 
uct Zg+i . . . Z„, we can divide the reduced state of the q 
particles in two ensembles. For each of these ensembles 
a different GHZ proof applies. If g = n, then (|8| must 
equal 1. Since these conditions cannot be fulfilled if X 
and Z of three or more particles have predefined values 
either —1 or 1, then we will conclude that the only hid- 
den instructions allowed in the LHV model are those in 
which X and Z of three different particles have not all of 
them predefined values. 

The challenge is to obtain the maximum possible detec- 
tion efficiency that can be reproduced with LHV models 
which satisfy (i) and (ii). Each of these LHV models is 
defined on a probability space (A,p), and is made up of 
subsets of instructions Ik,i.m C A, each of them charac- 
terized by three numbers: k is the number of particles for 
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which both observables {X and Z) are predefined (i.e., 
would give a detection when the observable is measured) , 
I is the number of particles for which only one of the ob- 
servables (X or Z) is predefined, and m — n — I — k is 
the number of particles for which none of the observables 
are predefined. 

Five lemmas are needed to prove our main result. 

Lemma 1. — In order to find the maximum detection ef- 
ficiency that can be reproduced with LHV models which 
satisfy (i) , it suffices to consider LHV models where each 
of the subsets Ik,i,m satisfies (i). 

Proof. — Suppose we find a LHV model compatible with 
a detection efficiency ry, such that some of the subsets 
Ik,i,m do not satisfy (i). Since the model must satisfy 
(i), we can always symmetrize it in all possible ways (by 
changing Z's to X's, and interchanging the different par- 
ticles) and consider an average of all these rearrange- 
ments. Clearly, this new model will have the same 77 and 
each of the Ik,i,m will satisfy (i). The new model will 
satisfy (ii) if and only if the original model satisfied (ii). 
■ 

Therefore, from now on we will only consider models 
such that each of the Ik,i.m satisfies (i). Each subset 
occurs with probcibility ^ Pk.i,7n 

In order to satisfy (ii) , the only subsets of instructions 
Ik.i,m allowed in our LHV models are those with k — 
0, 1, 2. In addition, the predefined values must satisfy (|8|, 
and the —1 and 1 values must be suitably distributed in 
order to refiect the fact that for the GHZ state all of the 
one qubit reduced density matrices are maximally mixed. 
Notice that these last two conditions are not particularly 
restrictive and can be easily satisfied. Therefore, in order 
to improve the clarity of the presentation, we will not 
insist on them hereafter. 

An upper bound on rj will follow from probabilistic 
considerations on each of the Ik,i.m- We will use the 
notation P/^ , ^ to refer to the probabilities of detection 
of the different variables within the sets Ik.i,m- 

Lemma 2. — The value of P/^ ^ ^ [Xi \X2, ■ ■ ■ , Xn) 
(and all the possible substitutions of Xi by Zi and per- 
mutations of the indexes) is n/ (2n — 2). 

Proof. — By definition. 



-f-r2.T.-2,o(^l|^2, ■ ■ • , Xn) 



-f/2,„-2,o(^2, • • ■ , Xn) 



(9) 



In the subset l2,n-2,o, only (2) instructions have prede- 
fined values for all the AT^'s. Since the total number of 
different instructions in /2,ri-2,o is (2)2""^, then 



Pi. 



, q{Xi, . . . ,Xn) 



1 



2«- 



(10) 



In order to calculate P/2,„_2,o(^2, ■ ■ ■ , Xn), we consider 
the subset S C /2,n-2,o where both Xi and Zi have pre- 
defined values. We also consider the complementary sub- 

and 



set - /2,„-2,o \ S. Clearly, P/,,„_,,,(5) 



Pl2,n.-2,o{^'^) 

we see that 



(s) ■ 



Reasoning in S and S'^ as above. 



_ _„(a:2,.. 

= Pl2.-n-2.a{^2, 
+ Pl2,^-2J,i^2, ■ 

1 C7^) 



, x„ 



, ^n\S)Pl2 „_2,oi^) 

.Xn\snPi2.,.^2Asi 

1 2n-2 



(2) 



(11) 



Lemma 3. — For every Ik,i,m different than /2,n-2,0) 
Pik I m(A^i|Ar2, . ■ . ,Xn) is either undefined or less than 
n/{2n - 2). 

Proof— li m > 1, P/^ ,^(Xi|X2, . . . ,X„) is not de- 
fined. If TO = 1, Pji^ , ^{Xi\X2, ■ ■ ■ ,Xn) = 0; hence, we 
consider only the case to = 0. In this case, there are only 
two subsets to consider, /i.,i_i,o and /o.n.o- 

Reasoning as in the proof of Lemma 2, 



-Pfo,n,o(^l|^2, ■ • • ,Xn) 



1 



and 



PIl,r^-l.oi^l\^2, • ■ • , Xn) 



2n- 1' 



(12) 



(13) 



which is always less than 2n-2 ■ ' 
Lemma 4- — For efficiencies higher than 77 = n/(2n — 2), 
there are no LHV models which simultaneously repro- 
duce all the quantum predictions ([s]) for q odd and 
S < q < n. 

Proof. — The LHV model must satisfy rj — P{Xi) — 
F{Xi\X2, Xn). The value of P(A:i|A:2, . . . , X„) must 
be less than or equal to the maximum of the values, where 
defined, of P/^ , ^(A'i|A'2, . . . , Ar„) for the different sub- 
sets Ik.i,m of the LHV model. According to Lemmas 2 
and 3, all of these values are less than n/(2n — 2). ■ 

Lemma 5. — For rj — n/(2n — 2), there are LHV models 
which simultaneously reproduce all the quantum predic- 
tions (jsj) for q odd and 3 < 9 < 

Proof. — We prove it by constructing explicit LHV 
models Mn{r]) = {{pk,i,m, Ik,i,m)} reproducing the quan- 
tum predictions for a given n and rj. Exact LHV models 
for n = 3, 4, 5 for rj = n/{2n — 2) are the following: 
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(!) = {(i,/2,i.o) , (^,/iA2) , (el. W)} , (14) 

^4 (I) = {(if ,^2,2,o) , (|:,/2,0,2) , (^,/l,0,3) , (^,^,0,4)}, (15) 

M5 (I) — {( ^2°°° '-^2,3,0) , ("^^i -^2,1,2) , ("^^^j -^2,0,3) , (^^,-^1,0,4) , (§11,-^0,0,5)} , (16) 



where, e.g., (||, /2,i,o) nieans that the model has instruc- 
tions /2,i,o with probability ||, etc. For higher n, we have 
obtained LHV models for rj = n/ (2n — 2) numerically for 
up to n = 15 qubits. For a given n and rj — n/{2n — 2), 
the LHV models are not unique. 

In addition, we have calculated numerically the maxi- 
mum background noise [5J [T^] as a function of the mini- 
mum detection efficiency required to violate the Mermin 
inequality. The results, for up to 7i = 8 qubits, are sum- 
marized in Fig. [l] 



the 5-partite Mermin inequality requires 5 spacelike sep- 
arated regions. When n tends to infinity, ?7crit tends to 
1/2, reflecting the fact that /o,„.o is a trivial LHV model 
compatible with the quantum predictions if < 1/2. 
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FIG. 1: Maximum (background) noise 1 — V as a. function 
of the minimum detection efficiency r) required to violate the 
Mermin inequality, when the state is 1/|GHZ„){GHZ„| -I- (1 — 
y)l/2", for n = 3,4, ... ,8 qubits. Note that if 1/ = 1, then 
ri = n/{2n - 2), and if ?? = 1, then V = 2'-^-"^/^. 

We have proven that a loophole-free Bell experiment 
using an n-qubit GHZ state and the correlations appear- 
ing in the n-partite Mermin inequality requires a detec- 
tion efficiency higher than n/{2n— 2). This result solves 
a long-standing open problem and is specially relevant 
for the 4 m], 5 Hg, and 6-qubit GHZ states [23] pre- 
pared in recent experiments. n/(2n — 2) is the threshold 
efficiency beyond which there is no LHV model which 
simultaneously satisfies all the quantum predictions Q 
and is the critical efficiency beyond which there is no 
LHV model reproducing all the quantum predictions for 
all the Bell inequalities, with two settings for q observers 
(3 < g < n) and one setting for the other n~q observers, 
contained in the Mermin inequality. This observation is 
of practical interest since, e.g., testing each of the (5) 2- 
2-2-1-1-setting Bell inequalities on a 5-qubit GHZ state 
requires only 3 spacelike separated regions, while testing 
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